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The fluctuations in an electron-ion plasma in thermal equilibrium are first treated for the collision- 
less case. Separate consideration is given to the contributions of the transverse and longitudinal 
modes, and two approaches are used; one in terms of the energy density of the thermally excited plasma, 
the other in terms of the permittivity e(/C,o>) appropriate to each mode. Some of the results have com- 
plete thermodynamic generality while some are based on the statistics of the plasma. Simple therm- 
odynamic arguments are used to derive the radiation resistance of short electric and magnetic dipoles 
in a collisionless plasma. The relation between the spatial correlation properties of the longitudinal 
mode and Debye-Hiickel theory is indicated. It is shown that in a plasma with collisions the fluctuation 
temperature is equal to the common electron-ion temperature for any degree of degeneracy and any 
velocity dependence of relaxation time. 



1. Introduction 

In an earlier paper [Burgess, 1941] 2 the author 
derived the electromagnetic field fluctuations of cavity 
radiation in a lossless nondispersive medium from 
the blackbody radiation energy density, and showed 
that an electric dipole immersed in the radiation would 
obey Nyquist's theorem. 

These results are now extended to a more general 
medium, specifically a plasma, which has a complex 
and dispersive permittivity and can support both 
transverse and longitudinal waves. The transverse 
mode, characterized by E and H being perpendicular 
to the wave vector K, was the only one involved in the 
earlier treatment. The longitudinal mode in which 
E and K are parallel and H is absent is accompanied 
by space charge density waves. In both modes damp- 
ing occurs, even in the absence of collisions, by virtue 
of the Landau mechanism. 

Since the dispersion equations are different for the 
two modes they will be considered separately; the 
thermal equilibrium energy density is calculated from 
the frequency-dependent permittivity and resolved 
into components arising from the fields and from the 
kinetic and potential energies of the particles. 

We are concerned throughout with an unbounded 
equilibrium plasma which is therefore spatially homo- 
geneous and statistically stationary. It is also assumed 



to be nonrelativistic; i.e., its temperature T is small- 
compared with mc 2 //c. Usually the ionic contribution 
to the fluctuations can be ignored if the ion mass is 
large compared with the electron mass but in the low- 
frequency spectrum for the longitudinal mode the 
ionic and electronic terms become equal. 

The spectral decomposition of a disturbance in the 
plasma is expressed in terms of the angular frequency 
o) and the wave number K; the linear response function 
of the plasma is then specified by the permittivity 
e(K, a>) which has different forms for the transverse 
and longitudinal modes. The presence of K in the 
permittivity represents the effect of spatial dispersion; 
an excellent discussion of the electrodynamics of 
plasmas involving both frequency dispersion and spa- 
tial dispersion is provided by the monograph of Silin 
and Rukhadze [1961]. 3 

Fluctuation spectra are derived using the stochastic 
current density given by the fluctuation-dissipation 
theorem applied to the imaginary part of the permittivity. 

Since in a vacuum cavity the electromagnetic field 
fluctuations are not zero it is necessary to demonstrate 
that the permittivity approach yields the correct result 
in the limit of Im (e) = 0. 

For a system in equilibrium, any temperature which 
can be used to characterize the properties must equal 
the plasma temperature; thus, for instance the temper- 
atures attributable to the electrons, the ions, and the 
electromagnetic field are all equal. Nevertheless 
erroneous statements have appeared in the literature 



'The research for this paper was supported by the Defence Research Hoard of Canada 
under Grant No. 9512-22. 

2 Burgess, R. E. (1941), Noise in receiving aerial systems, Proc. Phys. Soc. 53, 293-304. 
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to the effect that the fluctuation temperature of a 
plasma is equal to the electron temperature only in 
special circumstances e.g., a Maxwellian velocity 
distribution or a collision relaxation time independent 
of velocity. However in fact the equilibrium dis- 
tribution (Fermi-Dirac) for the electrons and ions suf- 
fices to ensure identity of these temperatures inde- 
pendently of degeneracy or of any velocity dependence 
of collision frequency. 

The temperature T enters the expressions derived 
in this paper through the fluctuation-dissipation 
theorem. Thus at frequencies small compared with 
kT/h the factor kT appears while at higher frequencies 
the more general form fico (exp Haj/kT— 1) _1 is appro- 
priate. For simplicity the kT factor will be used 
except where it is desired to integrate over all fre- 
quencies a spectrum which would otherwise yield a 
divergent result. In deriving the radiation resistance 
of a dipole immersed in a plasma it is of course im- 
material which form is used since the factor cancels 
in the comparison of mutually consistent expressions 
of Nyquist's relation and of the energy density. 

Throughout we shall use rationalized MKS units. 
Frequency spectra will be defined with respect to o> 
and exist for positive values; thus, for a fluctuating 
quantity X 



r 

Jo 



where co^ = 7i e 2 /me . This familiar result is a good 

approximation for wave numbers K small compared 
with the reciprocal of the Debye length D = 
(2n e 2 /eokT) 112 . In this case the Landau damping 
becomes very small and we write 

K/oj = O^oe) 1 ' 2 = (1 - a) 2 p l<o 2 yi 2 /c 
dk/dco = € (/We) 1/2 = (1 - co 2 Jw 2 )-^ 2 /c. 



The mean square electric and magnetic fields are 
found from the relation 



2 eE l = 2^ oH ^ 



(2.4) 



X 2 da)=X 2 . 



derived from Maxwell's equations and from the electric 
and magnetic energy densities 

while the kinetic energy density is 

^ = ^om(e 2 ^/m 2 co 2 ) = |(co> 2 )6o^;. 
The total energy density is thus 



2.1 



2. Transverse Mode 

Energy-Density Approach for the Transverse 
Mode 



In the transverse mode the electric and magnetic 
fields are perpendicular to the wave vector and no 
alternating component of space charge density occurs. 

The number of modes per unit volume in the wave 
number interval dK is 



dM = (kTK 2 /7T 2 )dK 



(2.1) 



taking into account two independent states of polariza- 
tion. The mean energy density thus has the spectrum 



U = (kTK 2 /7r 2 )(dK/da)). 



(2.2) 



It will here be assumed that the damping is suffi- 
ciently small for e and K to be substantially real for 
real co. It should be remembered that except for a 
few special functional forms of Re e(co) the imaginary 
part of e(co) would not necessarily be negligible at all 
frequencies. The requirement of causality imposes 
restrictions on the behavior of e(co) in the complex 
co-plane; one aspect of these is the Kramers-Kronig 
relations which connect Re e and Im e. 

Consider now a collision-free plasma for which the 
transverse permittivity is given by 



e/e = 1 — co 2 1 co 2 



We note that the effective permittivity is d{(oe)lda) 
for the nonmagnetic energy, since the sum of the elec- 
tric and kinetic energy densities can be written as 



U* + VS*=iEl(c+wkld»)- 



Finally, therefore the field fluctuation spectra of 
the total field component perpendicular to a given 
wave vector K are 

Yl=kTi^ 2 ^ 2 ^ 2 \TT 2 = kTo) 2 (l - (olJco^/Tr^eo (2.5) 

~Hl=kTa> 2 ^ 2 € 3 l 2 l7r 2 = kTo) 2 (l-o) 2 Jo> 2 )^ 2 l7T 2 c 3 fi . (2.6) 

The following relation between the three energy 
components exists 

f/el _ {/mag + jjkin = Lj] 



The total energy density spectrum has the frequency 
dependence 



U (0 = (kTl7r 2 c 3 )(o 2 (l - ollco 2 ) 1 ' 2 



(2.7) 



(2.3) for w > (Op. 



382 



If it is desired to integrate over all frequencies, kT 
must be replaced by Ji(jL>{cx\yfici)lkT — l)" 1 giving 



U-- 



7r 2 c 3 J 0) p exp 



,2X1/2 



*2 



(1ia>lkT) - 1 



dco. 



As the plasma density is reduced to zero (a) p — > 0) 
this integral increases to its maximum value given by 
the vacuum energy density of thermal radiation: 

C/ vac =77 2 ^ 4 r 4 /15c 3 ^ 3 . 

The fluctuation fields given above correspond to the 
resultant values. Thus the components in any direc- 
tion have mean square values equal to one-third of the 
given expressions. The mean square emf induced in 
a short electric dipole of effective length / is thus 

e&=PkTa?iLf€ l l 2 l3n 2 . 

This must equal (2/7r)kTRt as required by Nyquist's 
theorem, where Rt is the radiation resistance associ- 
ated with the transverse mode. Thus 



R* = Pv*ii%H l l 2 lbir. 



(2.8) 



Similarly a short magnetic dipole formed by a 
planar loop of area A has a mean square induced emf. 

Jfc = \ («M ) 2 ~Hl = I K*A*M- 

The radiation resistance for the magnetic dipole 
associated with transverse waves is thus 



/Jmag == ^2 ft| 4 /i 5/2 c 3/2/ 6irg 



(2.9) 



For both the electric and magnetic dipoles the radi- 
ation resistance is reduced by an ambient plasma com- 
pared with vacuum, the factors being respectively 
(e/eo) 1 ' 2 and (e/e ) 3/2 . 

This approach may also be extended to determining 
the spatial correlation of the field components. The 
cosine of the phase difference of the waves at two 
points (1, 2) separated by distance r is averaged uni- 
formly over all directions of arrival and all polarizations 
and gives 



E(ri)E*(rdJEl =(sin Kr)/Kr 

for the spatial autocorrelation coefficient appropriate to 
waves of wave number K = g)(jjlo€) 112 . Thus 



E( ri )E*(r 2 )„= (kT(oW 2 e^l7T*) (sin Kr/Kr) 

= (kT(DfjL {) J7r 2 ) sin [co(fi {) e) ll2 r]/r. (2.10) 

This result will be compared later with that derived 
using the e(K, co) method which includes the specific 
form of the imaginary part corresponding to Landau 
damping in a collisionless plasma. 



2.2. Permittivity Approach for the Transverse Mode 

In this approach the permittivity is taken to be a 
function of wave number as well as frequency in order 
to recognize the effects of spatial dispersion. For 
transverse waves the e(K, co) so defined determines the 
dispersion equation by the relation [Silin and Ruk- 
hadze, 1961, eq (6.8)]. 



e(K,o))le = K 2 c 2 /a) 2 . 



(2.11) 



For a Maxwellian equilibrium plasma with equal den- 
sity of electrons and (singly charged) ions the explicit 
form of the permittivity in the absence of collisions is 
[Silin and Rukhadze, 1961, eq (14.3)]. 



e(K, ft,)/ €o =l- (a> 2 p /(o 2 ) [Ffo+iirifiqe-*] 

- (fl 2 /co 2 ) [F(Q) +iir^Qe-^] 
where 

q 2 = o 2 m/2kT K 2 Q 2 = o) 2 M/2kT K 2 
while 



(2.12) 



F( 



?) = 2Xer*f* 



e**dy. 



For small values of X, F{X) -2* 2 -4¥ 4 /3+. . . for 
large values of*, F(X) ~ 1 + l/2* 2 + 3/4X^ + . . . 

When we are considering results obtained from the 
explicit form of e{K, w) we will usually be dealing with 
the case of large q where the imaginary part (Landau 
damping) will be relatively small and to a close approxi- 
mation 

Ree(K, co)/e - 1 - (<o 2 /w 2 ) (1 + kT K 2 /«> 2 m) 

-(D 2 /^ 2 ) (l+kTK 2 la> 2 M). 

The dispersion equation then becomes on further 
approximation 



co 2 - co 2 4- n 2 + K 2 (c 2 + kTIm) 



(0 2 p +K 2 C 2 



in view of the assumptions of large ion/electron mass 
ratio M/m and of a nonrelativistic plasma. 

The fluctuation current density normal to a given 
wave vector K has a spectral intensity determined by 
the plasma conductivity 

J 2 ~ = 2(2kTl7r) (-ulme). 

From Maxwell's equations the fluctuation field spec- 
trun are 

E*~= (UTIttu>) Im [l/(e-K 2 /"Vo)] (2.13) 

H*7=(4kTK 2 lTra>*(x® Im [l/(e-X 2 /<o>o)]. (2.14) 



383 



These expressions refer to the resultant quantity 
perpendicular to the wave vector and so give the sum 
of two independent orthogonal components. The Kco 
spectra may be integrated over K (or co) to yield spectra 
in terms of co (or K) and so give the purely temporal 
(or spatial) correlation functions. Performing these 
integrations invokes various sum rules which are listed 
in appendix 1 and which may be derived by contour 
integration. First, to find J$ we integrate the conduc- 
tivity of the Maxwellian plasma 



where 



o-(£, w)=— (olme(K, a>) = 7r 1/2 e co- 1 (w^e-« : 



for which 



/; 



+ OSQe-*. 



*(K 9 a>) dco= (7reo/2) (a>£+ftf). ( 2 .15) 



Thence 



J 2 K = 2n e 2 (kTlm + kTIM). 



(2.16) 



This expression may be derived directly by con- 
sidering the diffusion of the electrons and ions in the 
classical plasma. 

Of relevance in considering the relaxation process 
in a collisionless plasma is the autocorrelation of the 
current density: 



jK(to)J K (to + t) 



= 2kTn<>e*[m- x exp {-K 2 t 2 kTj2m) 

+ M" 1 exp(-^ 2 A:r/2M)]. 

This exhibits a gaussian autocorrelation function 
for the electrons and ions with respective relaxation 
times (m/kTK 2 )* and (M/kTK 2 )* corresponding to 
the gaussian thermal velocities of the particles. 

The electric and magnetic field components are 
found using the results (A2) and (A3) for integration 
over frequency 

E*=2kTMK,«>)=2kTl€o (2.17) 

Hl=2kTI,M>. (2.18) 

Integration over wave number to obtain the fre- 
quency spectra can be readily carried through if we 
ignore the dependence of e(K, co) on K when e appears 
in the form (e — K 2 la) 2 iJLo). This procedure is seen to 
be justified at small K (large q) since then we are 
neglecting kT/m by comparison with c 2 . On the other 
hand at large K (small Q) the contributions to the real 
and imaginary parts of € vary as l/K and are negligible 
compared with the K 2 /co 2 ijlo term. 

Thus taking e as a function only of a) we find an inte- 
gration over K space using integral (A8) 



p-ia = (D(fjL ey 12 . 



This result may be compared with the expression 
(2.10) for a lossless plasma; it has in addition an attenu- 
ation factor arising from the imaginary part of € (damp- 
ing), while the (sin fir)lr factor due to the phase dif- 
ference at the two observing points corresponds to a 
wave number oj/jl^ 12 Re (e 1/2 ). 

It is equivalent to the form given by Silin and Ruk- 
hadze [1961] in the second term of their equation 
(9.41). 

When the separation r goes to zero the total electric 
field fluctuation spectrum 



Ei = kTco 2 fjifRe(e 1 l 2 )/7T 2 . 



(2.20) 



When this is compared with the previous result (2.5) 
for a loss-free plasma it is seen that the appropriate 
generalization to a medium with damping (Im e ¥" 0) is 
to replace e 1/2 by Re (e 1/2 ). 

The electric field spectra derived in this section had 
their origin in the current density J which is determined 
by Im(— e). However it is seen that as Im(— e) tends 
to zero the electric field spectrum and its spatial 
correlation do not vanish but tend to the value obtained 
by considering the electromagnetic energy density in 
a cavity filled with a lossless medium. 

3. Longitudinal Mode 

3.1. Energy-Density Approach for the Longitudinal 
Mode 

The mode density per unit volume for the longi- 
tudinal mode is one-half of the expression used for 
the transverse mode because of the occurrence of 
only the single direction of the electric vector parallel 
to the wave vector. Thus the energy density spec- 
trum is 



U (0 = (kTK 2 /27T 2 )(dK/d(o). 



(3.1) 



E„(ri)Ei(r 2 ) - (kTa)Hol7r 2 )e-<* r (sin Mir 



The electric field may be derived from a potential 
and there is no associated magnetic field. The energy 
density comprises the contributions of the electric 
field and the kinetic and potential energies of the 
electrons, and we ignore the ionic contribution as 
being negligible when the mass ratio M/m is large: 

Un = ^ eo£ 2 + - n mul + - n mu 2 d (n 2 Jn 2 ) 

where & 2 is the spectrum of the electron velocity wave 
motion, u\ is the mean square thermal velocity SkT/m, 
and ra 2 is the spectrum of the electron density wave 
about the mean value n . The last term gives the con- 
tribution of the compressional energy density of the 
electron gas; since the pressure variations are given by 

Poj = ynJiT 
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and the compressibility by 

X=l/yn kT 
the density of the compressional energy is 



which may typically be very large compared with 
unity. 

The radiation resistance of a short electric dipole, 
of effective length /, associated with the longitudinal 
mode is scaled by the same factor and thus 



2PlX = ^ykT(n 2 Jn ) 



Kn*= l H^ 2 ™i*r. 



(3.3) 



and for longitudinal adiabatic compression we put 
7 = 3, yielding the result given above. 

By considering a longitudinal wave of the form represented by permittivity e(K, co). The dispersion 
exp i((ot — Kr), we can write the continuity equation equation for this mode is 



3.2. Permittivity Approach for the Longitudinal Mode 
Again the properties of the longitudinal mode are 



— iKrioUoj + icorio = 
and the equation of motion for unit volume as 
iwmnoUa) = — n^eE^ + iKmu^n^ 



e(£, a>) = 0. 



(3.4) 



For a Maxwellian plasma without collisions the 
permittivity assumes the form [Silin and Rukhadze, 
1961, eq (14.2)] 



where the last term arises from the gradiant of the e(K, Q))le =l—(Q) 2 l(jD 2 )2q 2 [F(q) — l + iTT ll2 qe q2 ] 
alternating pressure ynJiT. 



Thi 



conJKno = tt« = ieEJmoj(l — ulK 2 /o) 2 ). 



The total energy density is then 

"2 



U m = I €oE% [ 1 + <o> 2 + ulK 2 )/(co 2 - u 2 K 2 ) 2 ] . 



Introducing now the dispersion equation for longi- 
tudinal waves in its first approximation, for wave 
numbers small compared with the reciprocal of the 
Debye radius: 



j 2 = oj' 2 + K 2 u 2 



yields the results 



U„ = (<0*l<0*)6oEZ 



The reciprocals of the phase and group velocities are 

KI(o = (l-(o 2 p lw 2 y/ 2 lae 
dKJdu = (l-(o 2 p la> 2 )- 1 l 2 lu d 
leading to the result 



El = kTa* 2 (l - <o 2 l(o 2 yi 2 l27T 2 e ul 



- (Ct 2 laj 2 )2Q 2 [F(Q) - 1 + iir^Qe-W] (3.5) 

where q, Q, and F are as defined previously for the 
transverse mode. The factors in front of the square 
bracket terms have the common value 

(a> 2 l<o 2 )2q 2 = (n' 2 l(o 2 )2Q 2 = n e 2 /kTe K 2 = l/2K 2 D 2 

where 

D = (kTe l2n e 2 yl 2 

is the Debye shielding radius. 

For large values of q{KD <S 1) the asymptotic form 
of the permittivity on the assumption of M < m is found 
to be 

€(£, a)/* ~ 1 - (a> 2 /a> 2 ) (1 4- 3/2q 2 ). 

The corresponding dispersion equation takes on 
the well-known form appropriate to small wave 
numbers: 

a> 2 = a) 2 + (3kTlm)K 2 . 

The single Cartesian component of current density 
parallel to the wave vector has the spectrum 

Jl 0> =(2kTlTr)(-a>Ime). 
Integrating over frequency gives 



(3.2) 



Jl = kTe («>l + ft*) = n^JcTIm + kT/M) 



(3.6) 



This is larger than the electric field spectrum asso- u . , . , ,r c . u 1 . . , . , ., . 

•_^_j _-.i. .__ _j_ 1 .l: r__^„ which is one-hall ot the value associated with the trans- 

verse mode because of the two orthogonal components 
associated with the latter. Naturally the stochastic 



ciated with the transverse mode by the factor 

}(«£/<*Xc*/i4) 



current density which excites the plasma will have 
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three equal orthogonal components which cannot 
a priori be associated with any particular mode. 

This current density is related by the wave rela- 
tions to the other fluctuating quantities of interest 
namely: E, the longitudinal electric field; (/>, the electric 
potential, p the space charge density. For a wave 
of the form exp i{o)t — Kr) the relations are 



E = J/ia)€ 



<f> = EHK 



F 2 = I 2 



Me| 



<J>L 



= J|JKW|e| 2 



p=-iKeoE PL=/L^^ 2 kl 2 

so enabling the spectra of E, <£>, and p to be related to 
that of J. 
The electric field spectrum 



E% =(2*r/ir©)Im(l/€) 



(3.7) 



can be immediately integrated over all frequencies 
using (Al) to give 

E%=kT[ll€(K,co)-Ue(K,0)]. 

Now e(K, oo) = €o and e(K, 0) = e (l + l/K 2 D 2 ) 
whence 



E 2 K = kTl6 (l+K: 2 D 2 ). 



(3.8) 



In the evaluation of the static (zero frequency) permit- 
tivity e(K, 0) it is seen from (3.5) at q = Q = that the 
ionic contribution is equal to the electronic contribu- 
tion, and that the masses do not enter. In most other 
situations the large ratio M/m results in the ionic con- 
tribution to the permittivity being negligible; this is 
the case when q and Q are large compared with unity, 
which implies a wave number small compared with 1/D. 
The coulomb interaction in the ion-electron plasma 
vitally affects the electric field fluctuations in the 
longitudinal mode and this manifests itself through the 
appearance of the Debye distance providing a refer- 
ence length for the wave number K and for the spatial 
correlation. Integration over K gives 



E(n)E(r 2 ) = (kTI47T6 D 2 r)e- r l D = (n e 2 l27re 2 r)e- r l D . (3.9) 

This shows the Debye screening effect providing a 
cutoff in the coulomb interaction appearing in the spa- 
tial correlation of the electric field. 

The effect is shown in a different fashion in the 
potential fluctuations: 

~^ (t) = {2kTlTTK 2 a>) Im (1/c) (3.10) 

44 = kT/e K 2 (l + K 2 D 2 ) (3.1 1) 

0(r,)*(r 2 ) = (Ar/47r€o) (l-e"^)/r. (3.12) 



Here the potential correlation tends to a 1/r behavior 
at distances large compared with D. The mean square 
potential fluctuation at any point is 

$ 2 ~ = kmTT€oD = (kT/e) 2 l87rn D\ (3.13) 

The space charge density fluctuations are given by 



P L = (2£7X 2 e2/™)Im(l/e) 



P\ 



-kTeoK 2 l(l+K 2 D 2 ) 



(3.14) 



(3.15) 



pinWJ = - (kTeol*nD*)e-n»lr (r > 0) . (3. 16) 

This expression for the spatial correlation of the space 
charge excludes the singularity at r = which corre- 
sponds to the coexistence of a charge at r% if one is lo- 
cated to r\. The correlation function is negative due to 
the character of the coulomb interaction. 

If n and TV are the local densities of electrons and 
ions we can write (3.16) in the form 



(TV! - 7i,) (N 2 - n 2 ) = pip2le 2 =-n e r l D l2TrrD 2 (r > 0). 

Since the contributions of the electron-electron, elec- 
tron-ion and ion-ion interactions are dependent only 
on charge and not on mass we have for the individual 
spatial correlation of the density fluctuations 



AiViAAr 2 = AniA/i 2 =- A/i 1 A7V 2 =-n e-^/87rri) 2 (r>0). 

(3.17) 
This result gives the pair correlation in an electron- 
ion plasma compatible with the Debye-Huckel theory 
in its usual linearized form. 

The mean coulomb interaction energy per unit vol- 
ume of the plasma can now be found from these pair 
correlations: 



£/cou.- 2 J o 



(pip 2 /47re r)47rr 2 dr = — n e 2 /4<iT€QD. 



The kinetic energy density of the electrons and ions 
has a mean value of 

U kin = 3kTn . 

The ratio of the magnitudes of the coulomb and 
kinetic energy densities is 



|tfcoul|/£/kin = l/24W) 3 . 



(3.18) 



For the plasma to behave like a perfect gas, this 
ratio must be small which requires that n D 3 be larger 
than unity. Since D is proportional to n]J 2 this im- 
poses an upper bound on the required electron den- 
sity. An equivalent condition is obtained by demand- 
ing that the RMS potential fluctuation shall be small 
compared with kT/e which from (3.13) also requires 
that RqD 3 be larger than unity. Furthermore it is 
needed for the validity of linearizing the Poisson- 
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Boltzmann equations in the Debye-Hiickel theory. The autocovariance of N v is determined by the mo- 
Summing up, therefore, we may say that the present mentum relaxation time: 
discussion of coulomb interaction in the longitudinal 
mode is applicable only if cov [N v (to)N P (t + 1)] = (var N v ) exp [ — v{v)t\. 



n < (*7V2e 2 ) 3 . 



The z;-class current-density /^-spectrum has the auto- 
covariance 



4. Equilibrium Plasma With Collisions 
4.1. Current Fluctuation Spectrum 

Consideration will now be given to an equilibrium 
plasma with electron collision frequency v which, for 
generality of presenting the comparison with the col- 
lisionless case, will be assumed to be an arbitrary 
function of speed. 

The permittivity (to which the ionic contribution is 
neglected) is found from the Boltzmann equation ap- 
proach to be 



e(co) = 60 + 



3m 



Jo 



_d_ 

dv 



ia){io)-\- v\ 



4s7rCfo(v)dv. (4.1) 



Here fo(v) is the electron state occupation function 
and C = 2m s lh s with the electron density given by the 
integrated occupancy over all phase space 

r oc 

n = I 4<TrCv 2 fo(v)dv. 
Jo 

Now by the fluctuation-dissipation theorem 



Jv(to)Jr(t» + t) K = (e 2 v 2 /3V) cov \NMN v (t + 1)] 

where the factor 1/3 arises from considering a single 
component of J. 

Integrating over velocities gives 



J(toV(to + t)K=j J o "47rC/o(l -fo)v 4 e-*dv. (4.5) 

Using the cosine transform to give the frequency spec- 
trum yields 

JL = f£ \ X ^T~, 4vrCyi,( 1 - foWdv. (4.6) 
3tt Jo V* t a)" 



Now for electrons in equilibrium the occupation 
function / is the Fermi-Dirac function 



/o=[exp(|mt^-/x)/A7 T 4-l]- 



(4.7) 



where jx is the Fermi energy or chemical potential. 
We then have the identity. 



Jl (o = (2kTl7r)(-ojlme) 

2e 2 kT 
?>irm 



(4.2) 



dv~ kT Ml h) 



(4.8) 



I $ v I j^-f-^2 ) 47rC/o(v)di> (4.3) so that the current density fluctuations may be written 



for any one component of the current density. Inte- 
gration over frequency gives the sum rule. 



Pk 



: noe 2 kT/m = kTeoa) 2 



(4.4) 



which applies for any damping mechanism with any 
velocity dependence of the collisional relaxation time. 

The connection between the Boltzmann equation 
expression for the imaginary part of the permittivity 
and the current fluctuations is now made without ap- 
pealing to the fluctuation-dissipation theorem (4.2) by 
using statistical microscopic arguments. At this 
point the temperature does not appear and in fact it 
can only manifest itself when the statement that f is 
a thermal equilibrium distribution is inserted. 

We start by noting that in a volume V the variance 
of the number of electrons N v in the velocity interval 
(v, v + dv) is 

var N v = 4arCf Q. ~fo)Vv 2 • dv 

using the occupation fluctuation for Fermi particles. 



2e 2 kT 



Jo v--\-U) z dv 

Jo dv \v 2 -\-(x)' 2 



2e 2 kT 

377771 



IrrCfodv. (4.9) 



Thus the statistical approach and the thermodynamic 
result (4.3) agree as they must. The fluctuation tem- 
perature of the plasma and the kinetic temperature 
of the electrons are therefore identical and this iden- 
tity has been demonstrated for any arbitrary dependence 
of v on v and for Fermi-Dirac statistics. One can go 
further and consider a charged boson gas and the 
same result (4.9) is attained; to do this the + 1 is 
changed to — 1 in (4.7) and (1 — / ) is replaced in (4.8) 
by(l+/o). 

Statements which have appeared in the literature 
that the fluctuation and electron temperatures agree 
only under such restrictive conditions as a Maxwellian 
distribution or a velocity-independent relaxation time 
are seen to be erroneous. 
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5. Appendix A. Sum Rules and Spectrum 
Integrals 



U ft) e(K,a>) 2[e(K, 



1 



») e(A,0)J 



Jo ft) 

f=° dco 
Jo w 3 



1 



IT 1 



e(K,w)-£ 2 /ft>Vo 2e(K,°°) 



77^0 



(A-l) 

(A-2) 

(A-3) 



e(K,u))-K 2 l(o 2 fi IK 2 
J cM~ a> Im e) = (77/2)e (ft)2 + 0|) (A-4) 



f*Ksm_Kr 
Jo K 2 +^ 2 



^=^< 



(r>0) (A-5) 



I 



sin A> 7rx tt\— e~ rA 
aK=- 



o K(K 2 +^ 2 ) 



2 ^ 2 



(r^O) 



(A-6) 



Jo 


K 2 +A 2 


K 


sin Kr 


(K 2 +A 2 )(K 2 + B 2 ) 


K 3 


sin Xr 



sin /£r Ir , 77 ,„ , 

dK = -— A 2 e~ rA (r>0) 



(A-7) 



-,4r ,,-flr 



rfg=-^ ™ ,. (r>0) 



2 fi 2 -^ 2 



(A-8) 



(X 2 +.4 2 )(K 2 + fl 2 ) 



(A-9) 
(Paper 69D3-474) 
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